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ABSTRACT 



We have used a proper combination of multiband high-resolution HST- 
WFPC2 and wide-field ground based observations to image the galactic globular 
cluster NGC 6266 (M62). The extensive photometric data set allows us to deter- 
mine the center of gravity and to construct the most extended radial profile ever 
published for this cluster including, for the first time, detailed star counts in the 
very inner region. The star density profile is well reproduced by a standard King 
model with an extended core (~ 19") and a modest value of the concentration 
parameter (c = 1.5), indicating that the cluster has not-yet experienced core 
collapse. 

The millisecond pulsar population (whose members are all in binary systems) 
and the X-ray emitting population (more than 50 sources within the cluster half 
mass radius) suggest that NGC 6266 is in a dynamical phase particularly active 
in generating binaries through dynamical encounters. UV observations of the 
central region have been used to probe the population of blue straggler stars, 
whose origin might be also affected by dynamical interactions. The comparison 
with other globular clusters observed with a similar strategy shows that the blue 
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straggler content in NGC 6266 is relatively low, suggesting that the formation 
channel that produces binary systems hosting neutron stars or white dwarfs is 
not effective in significantly increasing the blue straggler population. Moreover, 
an anticorrelation between millisecond pulsar content and blue straggler specific 
frequency in globular cluster seems emerging with increasing evidence. 

Subject headings: Globular clusters: individual (NGC 6266); stars: evolution - 
binaries: close - blue stragglers 



1. Introduction 

In the last decades it has become obvious that some stars observed in Globular Clusters 
(GCs) are not the products of the passive evolution of isolated stars. Stellar encounters 
can produce binary systems and affect the evolution of the new produced binaries as well 
as any binaries that might have existed since the formation of the cluster. Thus the stellar 
content of a GC is intimately linked with the dynamical history of the cluster. In fact, 
the Color-Magnitude Diagrams (CMDs) of the core regions of GCs recently observed using 
HST-WFPC2 have shown the presence of a variety of exotic stellar objects, whose formation 
and evolution may be strongly affected by dynamical interactions. These include Low-Mass 
X-ray Binaries (LMXBs), (e.g. Heinke et al. 2001, Edmonds et al. 2003), Cataclysmic 
Variables (CVs) (e.g. Cool et al. 1998), and Blue Straggler Stars (BSSs) (e.g. the cases 
of 47 Tuc, Ferraro et al 2004; M3, Ferraro et al. 1997; M80, Ferraro et al., 1999a; M30, 
Guhathakurta et al., 1998, Ferraro et al. 2003a, hereafter F03). 

BSS were first discovered by Sandage (1953) in the globular cluster M3, and manifest 
themselves in the CMD as an extension of the Main Sequence (MS) above the Turn-Off (TO), 
mimicking MS stars with larger initial masses. There are two viable proposed mechanisms 
for BSS generation: the first is mass exchange in a primordial binary system and the second 
is the merger of two stars induced by stellar interactions (either single or binary) in a dense 
stellar environment. The fact that BSSs have been found in all properly observed stellar 
systems make them a powerful tracer of dynamical history of the parent cluster. 



1 Based on observations with the NASA/ESA HST, obtained at the Space Telescope Science Institute, 
which is operated by AURA, Inc., under NASA contract NAS5-26555. Also based on WFI observations 
collected at the European Southern Observatory (ESO), La Silla, Chile, within the observing programmes 
62.L-0354 and 64.L-0439. 
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Millisecond Pulsars (MSPs) may also be the by-product of collisions in the dense en- 
vironment of a cluster. They form in binaries containing a neutron star (NS) which is 
eventually spun up through mass accretion from the evolving companion. Since they are 
point-like objects and (in many cases) extremely stable clocks, they are invaluable tools to 
investigate the binary evolution in a very dense stellar environment (e.g. Rappaport et al. 
2001). It is important to study the overall properties of GCs hosting MSPs and to compare 
the distributions of the objects which may be significantly affected by dynamical encounters. 
Hence we have undertaken a long-term project to observe in the optical band the GCs which 
are rich in known MSPs, with particular emphasis on the BSS population. After having ex- 
amined 47 Tucanae (Ferraro et al. 2001, 2004; Mapelli et al. 2004) and NGC 6752 (Ferraro 
et al. 2003b, Sabbi et al. 2004), we now turn our attention to NGC 6266. 

NGC 6266 (M62) is a high density (logp c ~ 5.34), highly reddened (E(B - V) ~ 0.47) 
GC (Jacoby 2002; Possenti et al. 2003, hereafter P03) and is one of the most massive 
(M v = -9.19, Harris 1996 2 ). It has been classified possible Post Core Collapse (PCC) 
GC in the compilation by Djorgovski & Meylan (1993, hereafter DM93). Six binary MSPs 
have been recently discovered in this cluster (D'Amico et al. 2001a, 2001b; Jacoby et al. 2002; 
P03). NGC 6266 ranks fourth of the GCs in wealth of MSPs, after Terzan 5, 47 Tucanae and 
M15. Surprisingly, all MSPs in NGC 6266 are in binary systems (P03). As discussed in P03, 
the absence of known isolated MSPs in NGC 6266 cannot simply be ascribed to a selection 
effect since, for a given spin period and flux density, an isolated MSP is easier to detect than 
a binary MSP. If such a lack of isolated pulsars is not a statistical fluctuation, it must be 
somehow related to their formation mechanism and to the dynamical state of the cluster. 
Moreover, recent deep Chandra X-ray observations of NGC 6266 have revealed a very large 
number of X-ray sources — 51 sources were detected within the cluster 1.'23 half mass radius 
(Pooley et al. 2003), indicating that an high number of cataclysmic (and/or interacting) 
binaries should be present. These observational facts may indicate that NGC 6266 is now in 
a dynamical state where the rate lZf OTm of formation (and of hardening) of binary systems 
containing a neutron star and/or a white dwarf is much larger than the rate of disruption 
7^disr of such systems. 

In order to investigate in more detail the dynamical properties of this cluster, we have 
used HST-WFPC2 high resolution images of the central core and ESO-WFI wide field obser- 
vations to sample the more external regions. The combination of these two data-sets allows 
us to derive an accurate star density profile over the entire cluster extension. § 2 describes 
the datasets and the reduction procedures, while § 3 presents the derived CMD. § 4 describes 



2 For all references to Harris (1996) we have used the updated data set at the web site 
http: / / www. physics, mcmaster . ca/ Globular . html 
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the determination of the cluster center of gravity and the star density profile whereas § 5 is 
focused on the determination of the new cluster parameters. § 6 is devoted to the study of 
the Blue Straggler population in the central regions of NGC 6266 and, finally we discuss the 
results in § 7. 

2. Observations and data analysis 

The photometric data used here consist of two sets, (i) — The high spatial resolution 
set includes a series of HST-WFPC2 images obtained in August 2000, through the F555W 
(V), F336W (17), and F255W (UV) filters as part of a long term project (GO-8709, PI: F. 
R. Ferraro) aimed at studying the central stellar populations in a sample of GCs. In this 
dataset the planetary camera (PC, with the highest resolution of ~ O'.'046/pixel) is roughly 
centered on the cluster center, while the WF cameras (at a lower resolution of ~ &!l /pixel) 
sample the surrounding outer regions. (%%) — The wide field set — secured at the 2.2m ESO- 
MPI telescope at ESO (La Silla) in July 2000 using the WFI, which has exceptional imaging 
capabilities by providing a mosaic of 8 CCD chips (each with a field of view of 8' x 16') for 
a global field of view of 33' x 34'. The cluster is roughly centered on chip #2 and observed 
through the B, V, I broad band filters. 

2.1. Photometry and Astrometry 

The raw ground-based WFI images were corrected for bias and flat field, by using 
standard IRAF 3 tools. The point spread function (PSF) fitting procedure was performed 
independently on each B, V, and I image, using DAOPHOT II (Stetson 1994). A final 
catalog listing the instrumental B, V, and I magnitudes for all the stars in each field has 
been obtained by cross-correlating the three catalogs. The WFI catalog was finally calibrated 
by using the data-set by Rosenberg et al. (2000) for I magnitudes and by Brocato et al. 
(1996) for V and B magnitudes. 

The photometric reduction of the high resolution HST images was carried out using 
ROMAFOT (Buonanno et al. 1983), a package developed to perform accurate photometry 
in crowded fields and specifically optimized to handle under-sampled PSFs (Buonanno & 



3 IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Associa- 
tion of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 
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Iannicola 1989), as in the case of the HST-WF cameras. PSF-fitting instrumental magnitudes 
have been obtained using the standard procedure described in Ferraro et al. (1997, 2001). 
The final catalog of the F336W, F555W and F255W magnitudes was calibrated by using 
the zero-points listed by Holtzman et al. (1995). The HST F555W band photometry has 
been then transformed into the Johnson V system by using the stars in common between 
the WFPC2 and WFI catalogs. The photometric error of the final catalogs is dominated by 
the zero-point calibration uncertainties of ±0.05 mag. 

The Guide Star Catalog (GSCII) was used to search for astrometric standards in the 
entire WFI image field of view. Several hundred astrometric GSCII reference stars were 
found in each chip, allowing an accurate absolute positioning of the sources. An astrometric 
solution has been obtained for each of the 8 WFI chips independently, by using suitable 
catalog matching and cross-correlation tools developed at the Bologna Observatory . At the 
end of the entire procedure, the rms residuals are ps 0'.'2 both in RA and Dec. 

The small field of view (2'5 on the side) of the high resolution HST-WFPC2 images is 
entirely contained within the WFI chip #2, which imaged the cluster core. We used more 
than 2500 bright stars in the WFI catalog lying in the WFPC2 field of view as secondary 
astrometric standards, in order to properly find an astrometric solution for this catalog as 
well. We estimate that the global uncertainty in the astrometric solution is < 0'/3 both in 
RA and Dec. 

A master, homogeneous catalog of magnitudes and absolute coordinates including all 
the stars in the HST and the WFI catalog was finally produced. 

3. The CMD of NGC 6266 

Figure 1 show the CMD in the (V, U — V) plane for the entire HST sample. As can be 
seen all the cluster sequences are clearly defined and well populated. Particularly notable is 
the Horizontal Branch (HB) morphology. The CMD clearly shows an extended HB tail: this 
feature was already suspected on the basis of previous photometry (see for example Caloi et 
al. 1987, Brocato et al. 1996, and Contreras et al. 2005) but its true extension is revealed 
for the first time by the diagram shown in Figure 1 (see also Piotto et al., 2002). The HB tail 
extends ~ 1 mag below the MSTO and it is not uniformly populated. Similar morphology 
has been already observed in other clusters (NGC 2808 — Sosin et al 1997; M13 — Ferraro et 
al . 1997; M80— Ferraro et al. 1998; NGC 6752— Ferraro et al. 2003b, see also Catelan et 
al. 1998, Piotto et al. 1999 and Ferraro et al 1998 for a discussion on the reality of gaps 
along the HB). We defer the detailed discussion of the stellar distribution along the HB to 
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a future paper, where the HB morphology of these clusters will be compared to each other. 

NGC 6266 is in the sample of 62 GCs for which Ferraro et al. (1999b, hereafter F99) 
derived distance modulus and photometric properties. In doing this, F99 used the photome- 
try presented by Brocato et al (1996). As can be seen from their Figure 12b the HB appears 
quite dispersed when compared with the CMD shown here in Figure 1. Hence, since the zero 
point of the photometry presented here is based on the photometry by Brocato et al (1996), 
we used our new photometry in order to derive a more accurate level of the ZAHB (Vzahb)- 
Following the procedure described in F99, we obtained Vzahb — 16.25 ± 0.10, significantly 
(5V = 0.15) brighter than the values in F99. By adopting this new value, the reddening 
from Harris (1996) (E(B — V) — 0.47) and the metallicity from Carretta& Gratton (1997, 
hereafter CG97) ([Fe/H] CG97 = -1.07), we derived (from equation 4 in F99) My AHB = 0.68 
, an apparent distance modulus of (m — M) v = 15.57 ±0.15 and a true distance modulus of 
(m — M) = 14.11 ± 0.15 which corresponds to a distance of 6.6 Kpc. 

Recently Contreras et al. 2005 have found a large (> 200) population of RR Lyrae 
in NGC 6266. In order to determine the number of RR Lyrae variables lying in the field 
sampled by the HST-WFPC2, we cross-correlated our catalog with the list of variable stars 
by Contreras et al (2005). We identified 52 RR Lyrae in the HST/WFPC2 field of view. 
The positions of these stars in the CMD of Figure 1 are marked by large filled triangles. We 
also identified 68 RR Lyrae in the WFI catalogue. The mean magnitude of all 126 identified 
variables is < Vrrly >= 16.16 ± 0.4. By adopting the relation: 

M y (HB) = 0.22([Fe/H] + 1.5) + 0.56 

by Gratton et al (2003) and the metallicity of NGC 6266 ([Fe/H] CG97 = -1.07) we found 
My B = 0.65 and an apparent distance modulus (m — M)y = 15.55 ± 0.15, which is fully 
consistent with the value reported above. Hence in the following we adopt (m — M)y = 
15.57 ± 0.15, E(B - V) = 0.47 and (m - M) = 14.11 ± 0.15. 

4. The surface brightness and star density profile 

4.1. The center of gravity of NGC 6266 

For the first time star counts in the central region of the cluster can be computed using 
the CMD presented in Figure 1. Ferraro et al. (2003b, 2004) have shown that the knowledge 
of the position of individual stars in the innermost region of the cluster allows a high precision 
determination of the position of the center of gravity. To do this we applied the procedure 
described in Montegriffo et al. (1995) (see also Calzetti et al. 1993) and we computed 
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Cgrav by simply averaging the a and 5 coordinates of stars lying in the PC camera of the 
HST catalog. In order to evaluate any possible spurious effect due to incompleteness of the 
sample in the very inner region of the cluster, we considered samples with different limiting 
magnitudes (mpsssw < 20.5, 20 respectively). For each sample we computed the barycenter 
of the stars by using an iterative procedure (Ferraro et al 2003b). Both the determinations 
agreed within less than ~ 1". The center of gravity of the cluster (C grav ) turns out to be 
located at a J20 oo = 17 h 01 m 1278, 5 J200 o = -30° 06'46'.'0 with a typical la uncertainty of 075 
in both aj 2 ooo an d Sj 2 ooo, corresponding to about 10 pixels in the PC image. 

Figure 2 shows the 10" x 10" computer map around the cluster C grav as determined 
in this work (marked with a large cross at (0,0)). The cluster luminosity center derived 
by DM93 (small cross in Figure 2) is at about 3" NE from our determination, whereas the 
center of Harris (1996) is about 3" south. 

4.2. The star density and surface brightness profile 

The extended data set collected for this cluster offered the possibility to compute the 
radial star density profile. As first step, we define the star sample, paying particular care to 
avoid spurious effects due to possible incompleteness. Since such effects could be particularly 
important in the WFI sample, we restricted it to stars with r > 140". In the HST sample 
we have excluded the external region of the sample since the particular geometry of the 
WFPC2 field of view prevents an appropriate sampling of concentric annuli used to compute 
the radial density profile. A limiting magnitude of V ~ 20.5 was adopted for both samples. 
To summarize: 

1. the WFPC2/HST sample: all stars detected in the WFPC2 catalog with V < 20.5 and 
r < 93" from the cluster center; 

2. the WFI sample: all stars detected in the WFI catalog with V < 20.5 and r > 140" 
from the cluster center. 

Figure 3 shows the CMD in the (V, U — V) and (V, B — V) planes for the HST and 
WFI samples, respectively. As can be seen, the CMD of the WFI catalog (right panel in 
Figure 3) is severely contaminated by the disk population, which defines the almost vertical 
blue sequence (at 0.5 < B — V < 1.2) and, the bulge population at (B — V) ~ 1.5, which 
clearly indicates the presence of metal-rich stars. Figure 3 shows also a bright portion of the 
cluster Red Giant Branch (RGB) at V ~ 14 and B — V ~ 1.5 and, the blue extension of the 
HB at V ~ 16 and B - V < 0.5. 
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By using the combined data set shown in Figure 3 we computed the star density and 
surface brightness profiles, applying the standard procedure already described in previous 
papers (see Ferraro et al. 1999a, 2004). The entire photometric sample has been divided in 
41 concentric annuli centered on C grav , spanning a spatial range from 0" to 25'. Each annulus 
has been split in a number of sub-sectors (generally octants or quadrants, depending on the 
geometry of the field covered by the HST and WFI fields). The number of stars lying within 
each sub-sector was counted and star density was obtained by dividing the average star 
number by the area of the sub-sector. The surface brightness of each sub-sector has been 
computed by summing the luminosity of all the stars lying inside it, and normalizing the 
total luminosity to the subsector area. 

The average of the star density and surface brightness calculated in each subsectors, 
at a given radius, yields the stellar density and surface brightness of the parent annulus, 
respectively. The uncertainty in the average values for each annulus was estimated from the 
variance among the subsectors. Errors in the number counts can also be estimated from 
Poison counting statistics and are consistent with the empirically determined values. 

The star density values obtained for each annulus at different distance from the cluster 
center and the estimated errors are listed in Table 1 and shown in Figure 4. This is the most 
complete and extended density profile ever published for this cluster, since it samples the 
cluster population from the very inner core region out to r ~ 25' from the cluster center. 

The computed surface brightness can suffer relatively large fluctuations due to the small 
number statistics of the bright giants. Hence, we computed three radial profiles, removing 
the stars brighter than V = 12, 13, 14, respectively. Figure 5 shows that the overall structure 
of the profile does not change with the adopted magnitude limit. In contrast, the dispersion 
among different sub-sectors significantly decreases once the brightest stars are excluded. 

However, the surface profile shown in Figure 5 indicates a central brightness of fi v ^ ~ 
15.2 magarcsec -2 . This value is fully consistent with that listed by Harris (1996) (fi v ,o ~ 
15.35 magarcsec -2 ) and DM93 (/i„,o ~ 15.19 mag/arcsec 2 ), hence in the following we adopt 
this value for the central surface brightness of the cluster. 

4.3. The dynamical state of NGC 6266 

The shape of the radial density profile of GCs is generally used to infer the dynamical 
state of the system (see the pioneering work by King 1966 and Djorgovski & King 1986). In 
particular, Djorgovski & King (1986) defined two classes of clusters: the so-called King Model 
(KM) cluster and the PCC cluster, depending of the central density profile. KM clusters have 
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density profiles which can be reproduced by a single-component King Model (King 1966) 
with a flat isothermal core and a steep envelope characterized by two parameters: the core 
radius (r c ) and the tidal radius (r t ) or, alternatively, the concentration c = log(r t /r c ). The 
star density profile of a PCC cluster follows an almost pure power law with an exponent ~ 1. 
However Meylan & Heggie (1997) stated that the general dynamical status of a cluster can be 
deduced directly from the concentration value "c" , since all clusters with high concentration 
parameter value (c > 2) can be considered as collapsed or on the verge of collapsing. 

In order to model the observed profile, isotropic, single-mass King-models projected onto 
the cluster area have been computed by using the Sigurdsson & Phinney (1995) code. As is 
shown in Figure 4, the observed radial profile is well reproduced by a KM with an extended 
core (r c = 19"). Note that the value of the core radius is significantly larger (almost by a 
factor of 2) and the concentration (c = 1.5) significantly lower than r c = 10'.'7 and c = 1.7 
tabulated by Trager et al. (1995). Since the star density profile is not affected by statistical 
fluctuations due to the presence of a few bright giants, it offers the best route for deriving 
cluster parameters (see also Lugger, Cohn & Grindlay 1985 and Ferraro et al 2003b). 

Both the existence of an extended core and a relatively low value of the concentration 
parameter suggest that the cluster has not yet experienced the collapse of the core. Hence, 
the first result of this study is that NGC 6266 is not a PCC cluster. 

The star density profile shown in Figure 4 also allows us to evaluate the contribution 
of the background stellar population, which is substantial due to the fact that the cluster 
is located in the direction of the Galactic Bulge. Our star counts can be used to evaluate 
the projected Bulge star density down to V ~ 20.5. By averaging star counts in the region 
where the radial profile become flat (corresponding to r > 470" from the cluster center), we 
have estimated a value of ~ 92 stars arcmin" 2 . The background level is shown as a dashed 
horizontal line in Figure 4. The poorer fit of the King Model to the cluster in the region of 
transition to the background could arise from our modeling of the background as spatially 
uniform. 



5. The new cluster parameters 

Since the KM fitting of the observed stellar density profile of the cluster suggested 
values of r c and c significantly different from those previously measured, we used them to 
re-determine the cluster structural parameters. In doing this, we adopted the equations 
given by Djorgovski 1993 (from here D93). Following his assumptions the central luminosity 
density Pq is: 
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Po = s o/ (r c P) 

where r c is in parsecs and S is the central surface brightness in L VQ pc~ 2 , evaluated as 
follows: 

log(Eo) = 0.4 [26.362 - p y , (0)] 

and p is a function which depends on the cluster concentration (c), given by the following 
expression: 

log(p) = -0.603 x 10~ c + 0.302 

and /iy 5 o(0) is the central brightness of the cluster corrected for extinction. 

In computing these quantities we adopted the distance modulus and the reddening dis- 
cussed in Section 3. The resulting central luminosity density is log(pg) = 4.98 [L©pc~ 3 ]. 
If, according to D93, we use a mass-to-light ratio (M/L v ) = 3 for converting the cen- 
tral luminosity density into the central mass density, we derive a central density value 
of logpo — 5.46 [(M pc~ 3 )], which has to be compared with the value listed in D93 
logpo — 5.7 [(M Q pc~ 3 )]. All the derived parameters are listed in Table 2 along with the 
values by DM93 for comparison. 

To homogeneously compare the properties of NGC 6266 with those of other GCs, we also 
calculated new values of po f° r a sample of GCs previously studied with similar techniques 
(see F03, Sabbi et al. 2004, Ferraro et al. 2004). In doing this we adopted the distance 
scale defined by F99 (also the adopted reddening is from Table 2 by F99) and the r c derived 
from the radial density profile of each cluster (see F03, Mapelli et al. 2004, Ferraro et al. 
2003b). Results are listed in Table 3 along with the previous determinations by D93. The 
comparison shows that the different distance modulus adopted here (with respect to D93) 
has a direct effect on the derived value of the cluster central density, with a major impact 
on the derived binary-binary collision rate (see Section 7). 

6. The population of Blue Straggler Stars in the core of NGC 6266 

Special classes of stars, in particular those resulting from the evolution of binary systems, 
can be used as tracers of the GC dynamical evolution. This is the case of the BSS. From an 
observational point of view, the high angular resolution and UV imaging capabilities of HST 
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has allowed searches for BSS candidates in the core of highly concentrated GCs leading to 
a burst of activity in this sector (see for example Ferraro et al. 1997, 1999a, F03; Piotto et 
al. 2004). 

An obvious difficulty in studying BSS is the definition of an unequivocal selection cri- 
terion since the BSS sequence generally merges into the normal cluster MS without any 
discontinuity. However, previous work done by our group (Ferraro et al. 1997, 1999a, 2001, 
F03, 2004) has shown that the UV plane (the 771255, m 255 — m 336 plane in particular) is ideal 
for selecting BSS. In this plane the BSS sequence appears almost vertical and significantly 
brighter than the TO and Sub Giant Branch (SGB). In order to properly compare the BSS 
population in different clusters, here we follow the same procedure defined in Ferraro et al. 
1997 and adopted in F03. A limiting magnitude rriF255W — 19 has been adopted in selecting 
bright BSS once the CMD of each cluster has been shifted to match that of M3. Following 
this criterion, 47 BSS candidates have been selected in NGC 6266. Figure 6 (Left panel) 
shows the selection box adopted for the BSS in the UV plane. The position of the selected 
BSS in the (V, U — V^)-CMD is also shown in the Right panel of Figure 6. 

Since BSS are the result of two stars merging in a binary system, one can expect that in 
a dynamically relaxed stellar system, BSS are centrally concentrated with respect to normal 
stars (see Bailyn 1995) because of mass segregation. Indeed, the central relaxation time of 
even a sparse cluster is significantly lower than the lifetime of a 1-1.5 M MS star. Mass 
segregation has been observed in all the GCs properly observed up to now (see F03 and 
reference therein) except u Centauri (Ferraro et al. 2005), which may not be a true GC 
(Tsuchiya et al, 2004). 

When using UV-CMDs the obvious reference population is the HB, which appears quite 
luminous and well defined. HB stars have masses (~ 0.5 — 0.8 M & ) significantly smaller than 
those of BSS and their lifetimes are typically shorter than GC relaxation times. Figure 6 
(Left panel) shows the adopted selection box for HB stars in the UV plane. Indeed they 
are easely separable from the cooler (hence fainter) RGB stars. The right panel of Figure 6 
shows their position in the (V,U-V) CMD. As can be seen, the selection boxes defined in the 
UV plane assure the proper inclusion of the entire HB populations down to the faintest end 
of the branch. In the WFPC2 field of view, we have identified a total population of 395 HB 
stars (52 RRLyrae - see Section 3 and 343 non variable stars, respectively) 4 . 

We have determined the cumulative radial distribution of BSS and compared to HB stars 
as a function of the projected distance from the cluster center. Figure 7 shows the result. It 
is evident from the plot that BSS (solid line) are significantly more concentrated toward the 



Note that only non- variable HB stars have been plotted in both panels of Figure 6 
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center than HB stars. We applied a Kolmogorov-Smirnov test to both distributions to check 
the statistical significance of the detected difference: the probability that the BSS population 
and the HB stars are extracted from the same parent distribution is P = 2 x 10~ 4 . Thus, 
the central concentration of the BSS is confirmed with a > 99% statistical significance. 

6.1. BSS Specific Frequency 

In order to make a direct comparison with the BSS population in the previous published 
GCs (Ferraro 1997, F03, Sabbi et al 2004, etc), we computed the BSS specific frequency. As 
an example, Figure 8 shows the UV-CMD of NGC 6266 (left panel) compared to that of M80 
(right panel) presented in F03. From Ferraro et al. (1999a), the specific frequency (F^ s ) * s 
defined as the ratio between the number of BSS (Abss) to the number of HB (Ahb) stars. 
It's worth noticing that the number of stars in each post-MS phase is an excellent indicator 
of the sampled cluster luminosity (see Renzini & Buzzoni 1986). Hence F^ s i s a suitable 
tool for a quantitative comparison between BSS populations in different GCs. A further 
advantage is that F^ s is a purely observational quantity and can be easily computed in 
the UV-CMDs, where the HB population is well separated from the other sequences. The 
values of F|| s obtained in 6 GCs by F03, in NGC 6752 by Sabbi et al. 2004, in 47 Tucanae 
by Ferraro et al. 2004 and the one calculated in NGC 6266 are reported in Table 4, sorted 
according to increasing F^ s (the values of Abss an d Ahb are also reported) . All the clusters 
listed in Table 4 have been observed in the central regions by using the WFPC2/HST and 
hence represent an homogeneous sample for the sake of comparison. 

We found a particularly low BSS specific frequency {F^ s = 0.13) in NGC 6266, a 
value that lies between that of M13 and 47Tuc. It is interesting to notice that NGC 6266 
is the highest central density cluster in the sample but it has one of the lowest BSS specific 
frequencies. On the other hand F03 noticed that NGC 288, the cluster with the lowest 
central density in their sample, shows the largest BSS specific frequency. This is further 
evidence that the central density of a cluster is not the driving factor in determining the 
number of BSS (see also Davies et al., 2004). 

7. Discussion and Conclusions 

By adopting the cluster parameters computed in Section 5, we can now investigate some 
dynamical properties of NGC 6266. We computed the mean time interval between single- 
single (t ss ) and binary-binary (t bb ) collisions using equations 13 and 14 respectively from 



Leonard (1989): 



where a is the initial semimajor axis of each binary which can be estimated for each cluster 
according to the relation (F03) 

GM 

no is the central number density (number of stars per parsec 3 ) and V rms is the relative root- 
mean-square velocity, which can be approximated with the central velocity dispersion (er , 
see Table 5). Leonard (1989) derived his equation (14) under the assumption that the binary 
frequency is 100% in the core and the star mass is m* = O.5M . Assuming an average mass 
of = O.2M (Kroupa 2001) as suggested in F03, and = 0.5R & , we computed the 
number of single-single (ss) and binary-binary (bb) encounters occurring per Gyr (N SS;bb ) in 
the core of NGC 6266 to be: 

10 9 

N ss ,bb = [encounters Gyr 1 ]. 

tss,bb 

In Table 5 we give the results obtained for NGC 6266, compared to the 8 GCs which have 
been previously searched for BSS in the central regions by using the same observational 
strategy. The values of the encounter rate are normalized to those obtained for 47 Tuc. 
NGC 6266 has a predicted bb encounter rate N bb which is much greater than for any other 
analyzed GC. As shown in Table 5, the same result holds for the close encounters between 
two single stars leading to tidal capture. Both of these classes of dynamical interactions 
(together with the binary-single [bs] encounters) may lead to the merging of the two stars, 
thus potentially producing a BSS. Since all collision modes should be effective in increasing 
the population of BSS, one would expect to see a relatively larger specific frequency of BSS 
in NGC 6266. This expectation is not born out by the results of Table 4. The discrepancy 
might be reduced if the binary fraction in NGC 6266 were significantly smaller than in the 
other GCs, but this hypothesis can hardly be reconciled with the estimated high rate of 
binary formation 7Zf orm computed for this cluster (see later). Perhaps the results given in 
Tables 4 and 5 can be understood if the BSS in NGC 6266 result from primordial binaries 
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which have survived in the outer part of the cluster eventually wandering into the center 
where interactions drive them to merge (as we suggested for 47 Tuc, see Ferraro et al. 2004). 

On another side, the very high rate of binary-binary and binary-single star interactions 
in NGC 6266 may well explain why all of its known MSPs are in binary systems. In fact, P03 
suggested that NGC 6266 is in a particular dynamical phase in which the rate of binary sys- 
tem formation 7^f orm is larger than their destruction rate (IZdisr)- Following the prescription 
of Verbunt (2003), we recomputed those ratios for all clusters listed in Table 6. According 
to Verbunt (2003) 7?.f orm oc pl' 5 rl and 7^disr oc Po' 5 ?"" 1 , where po is the central density (ex- 
pressed in Mq pc~ 3 ) and r c is the core radius of the cluster. In Table 6 we give the values 
of 7Z{ OTm and TZdi ST normalized to 47 Tuc (as done by P03). A quite useful quantity is the 
ratio 7^form/7?-disr which quantifies the binary survival rate in each environment. NGC 6266 
shows, by far, the largest value of this ratio among the clusters listed in Table 6. The 2nd 
ranked cluster is M3 in which all the MSPs discovered so far are also in binaries. 

In summary, we have been able to recompute the ratio 7?.f orm /7?.disr and to show that 
(i) its high value nicely agrees with the hypothesis (supported by observations of MSPs and 
X-ray sources) that NGC 6266 has experienced (or is experiencing) a phase of very high 
production of binary systems. Additionally, we show that (ii) the highly effective binary 
production phase is occurring while the cluster has not yet undergone the collapse of the 
core, and that (in) the formation channel that produces a wealth of binary systems hosting 
neutron stars or white dwarfs is not efficient in producing BSS. 

We can also compare NGC 6266 with the sample of GCs searched for BSS with an 
approach similar to that used in this work. In particular, inspection of Tables 4 and 6 shows 
that in GCs hosting MSPs the specific frequency of BSS never exceeds 0.28. This sample 
includes the post-core collapsed globular cluster, NGC 6752, in which only a few BSS have 
been found (Sabbi et al. 2004). Conversely, up to now no MSPs have been found in M80 
and NGC 288, the two clusters with the largest specific frequency of BSS measured so far 
(probably generated by different channels: primordial binaries in NGC 288 and collisional 
binaries in M80). The difficulty in modeling the observational biases involved in discovering 
radio pulsars in GCs prevents any firm conclusion to date, but an anticorrelation between 
MSP content and BSS specific frequency seems emerging. GCs with large populations of 
BSS seem to have a small population of MSP. 
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Table 1. Surface density in different annuli around the 
NGC 6266 center. 



rl 


r2 


N 


err 


rl 


r2 




N 


err 


arcscc 2 


arcscc 2 


0.0 


2.0 


10.89 


0.04 


371.0 


404.0 
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026 


2.0 
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10.49 


0.04 
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437.0 





032 





025 
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845.0 





024 





004 


63.0 


73.0 


0.647 


0.025 


845.0 


899.0 





026 





003 


73.0 


83.0 


0.524 


0.012 


899.0 


953.0 





026 





021 


83.0 


93.0 


0.415 


0.021 


953.0 


1007.0 





025 





021 


140.0 


173.0 


0.107 


0.003 


1007.0 


1061.0 





025 





024 


173.0 


206.0 


0.081 


0.037 


1061.0 


1116.0 





027 





009 


206.0 


239.0 


0.061 


0.042 


1116.0 


1171.0 





027 





008 


239.0 


272.0 


0.045 


0.046 


1171.0 


1226.0 





029 





015 


272.0 


305.0 


0.038 


0.003 


1226.0 


1281.0 





029 





008 


305.0 


338.0 


0.037 
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Table 2. NGC 6266 parameters 





E(B-V) 


(m - M) 


^(0) 




r c 
[arcsec] 


c 


Log(po) 
[M pc- 3 ] 


DM93 


0.50 


13.70 


15.2 


13.64 


10.7 


1.8 


5.7 


This Paper 


0.47 


14.11 


15.2 


13.74 


19 


1.5 


5.46 



Note. — E(B-V) from Harris 1996 



Table 3. Structural paramenters for a sample of clusters searched for BSS 



Cluster [Fe/H] E(B-V) (M — m) Dist r c Log(p Q ) 

[Kpc] [arcsec] [M pc~ 3 ] 

This Paper DM93 



M 3 


-1.34 


0.01 


15.03 


10.1 


30 


4.0 


3.5 


M 13 


-1.39 


0.02 


14.43 


7.7 


40 


3.9 


3.4 


M 80 


-1.41 


0.18 


14.96 


9.8 


6.5 


5.4 


5.4 


M 10 


-1.41 


0.28 


13.38 


4.7 


40 


4.0 


3.8 


NGC 288 


-1.07 


0.03 


14.73 


8.8 


85 


2.3 


2.1 


M 92 


-2.16 


0.02 


14.78 


9.04 


14 


4.7 


4.4 


NGC 6266 


-1.07 


0.47 


14.11 


6.6 


19 


5.5 


5.7 


NGC 6752* 


-1.42 


0.04 


13.18 


4.3 


5.7 


5.6 


5.2 


NGC 104 


-0.70 


0.04 


13.32 


4.6 


21 


5.3 


5.1 



Note. - 



* Post core collapse cluster. 
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Table 4. BSS population in a sample of 
GGCs. 



Cluster 


Nbss 




pBSS 


1 1/2 I'c 


M13 


16 


228 


0.07 


1.15 


NGC 6266 


47 


395 


0.12 


0.63 


NGC 104f 


53 


314 


0.17 


0.87 


NGC 6752* 


28 


156 


0.18 


0.82 


M10 


22 


81 


0.27 


0.85 


M3 


72 


257 


0.28 


0.73 


M92 


53 


160 


0.33 


1.07 


M80 


129 


293 


0.44 


1.07 


NGC 288 


24 


26 


0.92 


0.71 



Note. — | In 47 Tuc the selection of BSS 
and HB has been performed using the CMD 
plane m F2 i 8 w ) m F2i8W ~ ?™F439W instead of 
the m F 255W,"^F255W - "^F336W plane. 

Note. — * Post core collapse cluster. 
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Table 5. Expected number of single star (N ss ) 
and binary-binary (N bb ) encounters per Gyr 



Cluster 


a hs (AU) 


oo 


N ss 


N bb 






[km s^ 1 ] 






M13 


0.39 


7.1 


0.07 


0.17 


NGC 6266 


0.10 


14.3 


3.19 


2.07 


NGC 104 


0.15 


11.5 


1 


1 


NGC 6752 


0.98 


4.5 


0.16 


1.04 


M10 


0.45 


6.6 


0.05 


0.11 


M3 


0.63 


5.6 


0.16 


0.70 


M92 


0.57 


5.9 


0.28 


1.04 


M80 


0.13 


12.4 


0.33 


0.29 


NGC 288 


2.36 


2.9 





0.02 



Table 6. MSP content, destruction and formation rate 



Cluster 


Isolated 
PSRs 


Binary 
PSRs 






7^-form/ 


M13 


2 


3 


6.51 x 10~ 2 


7.18 x 10~ 2 


1.1 


NGC 6266 





6 


2.65 


1.09 


2.97 


NGC 104 


7 


15 


1.00 


1.00 


1.00 


NGC 6752 


4 


1 


0.18 


6.72 


3.21 x 10~ 2 


M10 








4.61 x 10~ 2 


0.14 


0.40 


M3 





4t 


0.11 


8.83 x 10~ 2 


1.58 


M92 








0.21 


0.47 


0.56 


M80 








0.51 


1.97 


0.32 


NGC288 








1.6 x 10~ 3 


4.72 x 10~ 3 


0.41 



Note. - | For 3 of these pulsars no timing/orbital solution has 
been published yet. 
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Fig. 1. — CMDs for stars in the HST/WFPC2 FoV. RRLyrae stars lying in HST/WFPC2 field of view are 
marked with large filled triangles. 
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Fig. 2. — Computer map of the inner 10" x 10" region of the cluster. The large cross at (0, 0) coordinates 
indicates the adopted C grav . The Ci um by Djorgovski (1993) and Harris (1996) are labeled with D and H, 
respectively. 
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Fig. 3. — CMDs for stars in the two samples. Left panel: the high resolution HST catalog (r < 93")in 
the (F555W, F336W - F555W) plane. Right panel: the wide-field WFI catalog in the (V, B — V) plane, 
only stars with r > 140" from the cluster center are plotted. Stars with V > 20.5 (shaded region) have been 
excluded from the construction of the density profile. 
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Fig. 4. — Observed radial density profile with respect to the the adopted C grav . The solid line is the best fit 
King model (c = 1.5) to the observed density profile over the entire extension. A value of ~ 92 stars arcmin~ 2 
is representative field contamination (dashed horizontal line) . The adopted value for the parameter Wo (the 
central potential parameter as defined by King 1996) is also reported. 
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Fig. 5. — Radial brightness profiles computed, for the WFPC2/HST sample, by removing the stars brighter 
than V = 12, 13, 14, respectively, the faintest limit is V = 20.5 (see Figure 3). The solid line is the best-fit 
solution. 
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Fig. 6.— BSS (open circles) and HB (open squares) in the (F255W, F255W - F336W) CMD of NGC 6266 
(left panel) and (F555W, F336W - F555W) CMD (right panel). Note that only non variable HB stars are 
plotted. The adopted selection boxes in the UV-CMD are shown. 
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Fig. 7. — Cumulative radial distribution of BSS (solid line) with respect to the HB stars (dashed line) as 
a function of their projected distance (r) from the cluster center. The probability that the two populations 
are extracted from the same distribution is P = 2 x 10~ 4 
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Fig. 8.— BSS (open circles) in the (F255W, F255W - F336W) CMDs of NGC 6266 (left panel) and M80 
(right panel) for comparison. The CMDs have been shifted as suggested in F03. The horizontal line shows 
the limiting magnitude of the adopted selection. 



